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In this Comment I argue that Silveirinha's criticism of my earlier work [M.G.Silveirinha, Phys. 
Rev. B, 80, 235120 (2009)] is based on an unphysical excitation model which involves an "external 
current" which overlaps with a continuous medium but is not subject to constitutive relations. When 
this excitation model is replaced by the conventional and experimentally-relevant model of excitation 
by external electromagnetic fields, it can be easily shown that Silveirinha's formulas contain the very 
results he wanted to disprove. This and a few other misconceptions present in Silveirinha's paper 
are subject of this Comment. 



I. INTRODUCTION 

This paper contains a comment on the paper! by 
M.G. Silveirinha entitled "Poynting vector, heating rate, 
and stored energy in structured materials: A first- 
principles derivation" [Phys. Rev. B 80, 235120 (2009)]. 

As is customary for all APS journals, the Comment was 
first sent to Silveirinha for review. Silveirinha has written 
an 11-page report on my Comment and recommended 
its rejection. I disagreed, the Comment went through 
the ordinary process of peer-review and was eventually 
rejected. Still, I did not agree with the scientific argu- 
ments used to justify the rejection. I have then reworked 
the Comment into a stand-alone paper with independent 
conclusions. This paper— is now published and is entitled 
"On the current-driven model in the classical electrody- 
namics of continuous media" [J. Phys.-Cond. Matt. 22, 
485401 (2010)]. Although the above publication makes 
all the essential points, I still feel that the original Com- 
ment should be made available to the community as it 
is, in some respects, more direct and specific. 

Finally, it should be mentioned that a different com- 
ment^ on the same paper— by Silveirinha has been sub- 
mitted and published. 



II. THE ORIGINAL COMMENT 

Recently, Silveirinha has criticized^ several earlier pa- 
pers^—, two of which are authored by myself. In par- 
ticular, Silveirinha has stated that the conclusions made 
by me in Ref. [3, as well as similar conclusions made by 
Richter et al. in Ref. HI, "are founded on fundamental 
misconceptions and mistakes" . To support this claim, 
Silveirinha has offered a rather lengthy derivation which, 
according to Silveirinha, disproves the results of Refs. |4j 
and [5j. However, I have found that Silveirinha's deriva- 
tion itself is based on a few misconceptions. The chief 
among these is the contention that an electromagnetic 
wave with a purely real wave vector can propagate in 
absorbing materials (either "meta-" or otherwise), that 
such non-decaying waves can be obtained by imposing an 



"external current" which overlaps with the medium but is 
not subject to constitutive relations, and that solutions 
thus obtained can be used to derive information about 
the rate at which the medium is heated by radiation (the 
heating rate). This and a few other "misconceptions and 
mistakes" which are present in Silveirinha's work are sub- 
ject of this Comment. Gaussian system of units is used 
throughout. Equations of Silveirinha's paper are refer- 
enced as (S#). 

A) The assumption of purely real wave vec- 
tor and the external current. Most of Silveirinha's 
derivations, including all derivations that lead to results 
of any practical significance, are carried out for Bloch 
waves with a purely real wave vectors k. Silveirinha is, 
of course, aware that natural Bloch modes in periodic 
media with some amount of absorption are, necessarily, 
decaying waves. In a typical experiment, these decay- 
ing modes are excited by external radiation. The decay 
is mathematically represented by the imaginary part of 
k. Thus, if k is taken to be purely real, the medium is, 
by definition, non-absorbing. Calculation of the heating 
rate in such a medium is meaningless: under the con- 
dition Imk = 0, any reasonable calculation must yield 
zero. 

Silveirinha, however, contends in his paper that a real- 
valued Bloch wave vector k is not incompatible with 
losses. He argues that one can consider an external cur- 
rent rather than an external radiation as the source for 
the electromagnetic wave in the medium. This is quite 
fine as long as the current and the medium do not over- 
lap. But Silveirinha requires that the external current 
be of the form j e = Aexp(zk • r), where A is a constant 
amplitude, so that j e is nonzero in the region of space 
occupied by the medium. It must be emphasized that 
the current j e introduced by Silveirinha is not created by 
the charges associated with the medium, does not obey 
any constitutive relations (and, therefore, is independent 
of the fields existing in the medium), and, presumably, 
can be manipulated by the experimenter at will. 

The excitation model proposed by Silveirinha and de- 
scribed briefly above is unphysical because it can not be 
realized experimentally. The obvious question Silveirinha 
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should have answered is the following: how can an exter- 
nal current of the form j e = A exp(zk-r) be physically cre- 
ated in the medium? If we set aside heating, mechanical 
stress, acoustic waves and other similar influences (which, 
anyway, are highly unlikely to generate the required cur- 
rent), we are left with the possibility of applying external 
electromagnetic fields to the medium. But any current 
produced in such an experiment is subject to the usual 
constitutive relations and is, therefore, a part of the in- 
duced current which is denoted in Ref. [l| by j^. Another 
option is to insert multiple wires into the medium and 
to run externally-controlled currents through the wires. 
This will, of course, destroy the medium and no reason- 
able experimentalist would propose or try to implement 
such an excitation scheme. It is also highly unlikely that 
pre-determined spatially-varying currents can be main- 
tained in the wires by external voltages. 

Thus, the idea of exciting the medium with a pre- 
determined external current which overlaps with the 
medium is unphysical and can not be used to derive any 
useful physical quantity. But even from a purely formal 
point of view, any results derived from such an excita- 
tion model can not be used to criticize Ref. |j because 
in this paper I have clearly and repeatedly stated that I 
consider only the case when all currents in the medium 
are induced and there are no external or "free" currents 
within the medium. (Note that the conductivity current 
is included in the induced current in this paper and in 
Ref.0.) 

B) Heating rate with and without the assump- 
tion of a purely real wave vector. Now let us con- 
sider how Silveirinha computes the heating rate. The 
starting point will be Eq. (S60). To save a few notations, 
I will re-write this equation in Gaussian units and omit 
all subscripts. Then Eq. (S60) takes the form 



= — Im 

8tt 



E* • eE - 



E* • k x (/T 1 — 7)k x E 



(1) 

Here tensors (dyadics) are denoted by a hat, the symbol 
" •" is reserved for a dot product of two vectors, no spe- 
cial symbol is used to denote the action of a tensor on 
a vector (for example, eE), and all operations are eval- 
uated from right to left; E is the macroscopic electric 
field and I is the unit tensor. All fields are assumed to 
be monochromatic and the common exponential factor 
exp(— iu>t) is suppressed. Note that I have used the first 
identity in (S36) to transform (S60) to the form (p}. In 
fact, ((T|) is equivalent to (S60) but is written using more 
conventional tensor notations. 

Silveirinha evaluates expression ([T]) as follows. He uses 
the vector identity a • (b x c) — (a x b) • c (equivalent to 
the second identity in (S36)) to re- write the second term 
in the square brackets as 



(E* x k) • [(/T 1 - Z)k x E] 



Then Silveirinha uses Eq. (SlOa), namely, k x E = 
(u;/c)B. This equality follows from the macroscopic 
Maxwell's equations applied to a plane wave. It can be 
seen that the factor k x E in the end of expression @ can 
be replaced by (w/c)B. However, the expression E* x k 
in the beginning of this expression can only be replaced 
by — (oj/c)B* if k is purely real. Silveirinha makes this 
assumption about k and transforms © to the form 



b* • (/r 1 - j)b 



(3) 



The term proportional to I is then omitted since its imag- 
inary part evaluates to zero, the field B is expressed 
in terms of the field H using the constitutive relation 
B = jftH and Silveirinha arrives at his Eqs. (S61) and 
(S62). The important point here is that the transition 
from (S60) to (S61) and (S62) requires that Imk = 0. 

But if we set aside the unphysical model in which the 
medium is excited by an external current which overlaps 
with the medium, then k is not a mathematically inde- 
pendent variable but is determined from the Maxwell's 
equations and constitutive relations. The equality Imk = 
is only possible if Ime = Im/i = 0. Therefore, the 
final result for the heating rate derived by Silveirinha, 
Eq. (S62), evaluates to zero under the assumption that 
was used to derive it. It can be concluded that the result 
(S62) is not a "first-principle, completely general" deriva- 
tion of some useful physical quantity, as Silveirinha has 
claimed, but a lengthy and laborious proof of the identity 
= 0. 

It is possible, however, to evaluate the expression ([1]) 
differently without making any assumptions about k. To 
this end, we use the two Maxwell's equations 



k x E 



k x H = e"E 

c 



(4) 



Here I have assumed that all currents in the medium are 
induced and subject to constitutive relations. From ((4]), 
we can also obtain 



k x /r x k x E = - (-) 2 £ E 



(5) 



Substitute this expression into ([T]). The term propor- 
tional to e will cancel to yield 



q = - — (-) 2 i m (e* ■ k x k x E) 

87T \W/ 



(6) 



Next, use the identity k x k x E = k(k • E) - fc 2 E to 
obtain 



(2) 



q= 8^U Im [|E| 2 fc 2 -(E-k)(E*.k)] . (7) 

This expression is equivalent to the one derived by me 
earlier in Ref. H] (Eq. (54) of this reference) for the case 
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of general anisotropic nonlocal media. If the medium 
is isotropic, it can support only transverse waves whose 
wave number is k 2 — (uj/c) 2 e[i with e and fx being scalars. 
Then ([7]) is simplified to 



qjjE| 



■Im(e^) 



(8) 



Again, this result was derived by me in Ref. |3 (see 
Eq. (32)). 

Thus, it can be seen that Silveirinha's Eq. (S60) con- 
tains, in fact, the very results he wanted to disprove. The 
only reason Silveirinha has obtained a formula which ap- 
pears to be different from or © is because he has 
used a method to evaluate his equation (S60) which is 
only valid when Imk = 0. My method of evaluating 
(S60) makes no assumptions about k. It is applicable, in 
particular, when Imk = 0. In this case, ([7]) and (HJ, as 
well as Silveirinha's results (S61), (S62), all evaluate to 
zero and are, in this sense, equivalent. But, unlike Sil- 
veirinha's results, formulas ([7]) and ([5]) can also be used in 
the physically interesting case of a complex wave vector 
k. 

Silveirinha would, perhaps, argue that, for complex 
wave vectors, (S60) is itself invalid because its deriva- 
tion requires that the second term in the right-hand side 
of (S58) be neglected. This term was shown to be ex- 
actly zero when Imk = 0. However, in electromagnet- 
ically homogeneous media, this term is negligibly small 
anyway. This is evident from the following consideration. 
The tensor G p o defined in (S57) is real and Hermitian as 
long as k ■ (r — r') <C 1. Since the double integration 
in (S58) is over the volume of the same elementary cell, 
the maximum phase shift k • (r — r') is of the order of 
or smaller than \k\h, where h is the characteristic cell 
size. For the medium to be considered electromagnet- 
ically homogeneous, this phase shift must be negligibly 
small. Therefore, the result of double integration in (S58) 
is purely real with the same precision as the precision of 
the underlying approximation. In other words, the er- 
rors that a made when the original periodic medium is 
replaced by a spatially-uniform medium with effective pa- 
rameters are of the same order of magnitude as the errors 
which are made when the last term in (S58) is neglected. 

C) The lack of equivalence between magnetiza- 
tion and nonlocality. Another serious misconception 
perpetrated by Silveirinha is that, at high frequencies, 
magnetic response of matter can be understood and is 
physically indistinguishable from nonlocality in the di- 
electric response. To be fair, we should note that Sil- 
veirinha did not come up with this misconception on his 
own. The misconception can be traced to the highly re- 
spected book by Landau and Lifshitzi and was recently 
popularized by Agranovich 8 . However, it is elementary 
to see that no such equivalence exists. 

Indeed, consider the induced current in a medium with 
electric and magnetic polarization. Let, for simplicity, 
the medium be isotropic with scalar coefficients e and fi. 



For monochromatic fields, the induced current is 



-iwP + cV x M 



(9) 



where P = XeE and M = XmB, % e = (e — l)/47r and 
Xm = (m — l)/47r/i. Traditionally (e.g., in Landau and 
Lifshitz's book or in Agranovich's review), the expression 
© is evaluated for a plane wave with the wave vector k 
by writing 



V x M = V x XmB = Xm V x B = i-Xmk x k x E . (10) 

Thus, it is found that the current can be equivalently 
written as 



e n0 ni(w,k) - I 
> 

47T 



E 



(11) 



where 



U^k)=£l-(-) 2 ^kxkx . (12) 

Upon observing that exactly the same induced current is 
obtained in a medium with permeability equal to unity 
and nonlocal permittivity e n0 ni(w, k) as the current ob- 
tained in a medium with local permittivity e / 1 and 
local permeability fj, ^ 1, the conclusion is drawn about 
the physical indistinguishability of the effects of magne- 
tization and nonlocality of the dielectric response. 

The above line of reasoning is erroneous because V x 
XmB ^ XmV x B. The correct formula is 



V X XmB = XmV X B+ (VXm) x B . (13) 

The last term in the above equation is often lost sight 
of because it is assumed that Xm = const and, conse- 
quently, Vxm = 0. However, even if the medium is quite 
uniform spatially, it always has a boundary, and at the 
boundary the function Xm(r) experiences a jump. Dif- 
ferentiation of this jump gives rise to a surface current 
which can not be accounted for in the model of nonlo- 
cal permittivity (TT2"j) . To be more precise, accounting for 
the surface current within the spatial nonlocality model 
would require the use of position-dependent and highly 
singular integral kernels in the influence function which 
would contain derivatives of the delta function. This pos- 
sibility is never considered. In particular, in Refs.@H the 
surface current is simply ignored. 

Thus, the deficiency of the equivalence model is that 
it places the emphasis on the dispersion relation but ig- 
nores the medium boundary. Although infinite bound- 
less media can be considered as purely mathematical ab- 
stractions, such media do not exists in nature. And a 
dispersion relation alone does not characterize a finite- 
size sample: the latter also has an impedance. Indeed, 
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the dispersion relation which follows from (|ll[) . (TT2t is the 
familiar equation k 2 = (ui/c) 2 efi, where I have used the 
fact that waves in isotropic media are transverse. This 
dispersion relation is invariant with respect to the trans- 
formation e — > /3e, /i — > where /3 is an arbitrary 
complex constant. In particular, we can take /3 = fi. It 
then would appear that a medium with some values of e 
and is indistinguishable from a medium with s' = s/i 
and fjf = 1. But the two media have different impedances 
and, therefore, different transmission and reflection prop- 
erties. 

D) The meaning of homogenization. Yet another 
misconception encountered in Silveirinha's work is re- 
lated to the notion of homogenization. Silveirinha seems 
to believe that homogenization is equivalent to field av- 
eraging. However, fields can always be averaged but 
not every medium is electromagnetically homogeneous 
at a given frequency. The goals of homogenization are 
the following. Firstly, one needs to prove that solutions 
to the macroscopic Maxwell's equations in a nonmag- 
netic spatially-inhomogeneous medium characterized by 
the "true" permittivity etrue(r) converge in some phys- 
ically reasonable norm to solutions in a medium of the 
same overall shape but with spatially-uniform effective 
parameters e e ff and /z e ff. Secondly, the effective parame- 
ters must be computed explicitly in terms of £true( r )- The 
original and the effective media must be "electromag- 
netically similar" for varying incident fields and varying 
shapes of the sample. 

None of this has been done in Ref. 1!. Silveirinha did 
not consider the conditions under which spatial averaging 
of fields is physically meaningful, nor did he derive the 
effective parameters in terms of £true(r). True, more was 
done by Silveirinha in Ref. [^, where a numerical proce- 
dure to determine the effective parameters was proposed 
and implemented. However, this development was not 
used in Ref. [l|. Besides, Ref. HI is affected by the same un- 
physical excitation model as Ref. Q} the Bloch wave num- 
ber k is viewed in both references as a mathematically- 
independent variable "imposed" by the external current 
although it is well known that to solve a photonic crys- 
tal, k must be computed self-consistently. In any case, 
Silveirinha's derivation presented in Ref. .lj is simply a 
chain of definitions which, after some terms are neglected, 
result in familiar macroscopic formulas, such as (S60) 
(which Silveirinha then evaluates incorrectly). These for- 
mulas contain effective parameters which, at that point, 
are purely phenomenological and which have not been 
computed either "from first principles" or in any other 
meaningful way. For example, in Eq. (S60), there appear 
some effective medium parameters e r and /i r . How are 
these related to £tmc(r)? How does Silveirinha know that 
Hr is different from unity? 

In this respect, an interesting result has been recently 
reported by Menzel et al^ It was found that a typical 
"metamaterial" can not be reasonably characterized by 
effective parameters in the spectral region in which mag- 
netic resonances are excited. This result indicates that, 



contrary to Silveirinha's assertion, "artificial magnetism" 
can not be obtained in electromagnetically homogeneous 
materials. Here, following Ref. [Il| and a naive analogy 
with natural magnetism, Silveirinha has confused mag- 
netic polarizability of subwavelength particles with mag- 
netic permeability of a periodic medium made of such 
particles. In the case of natural magnetism, it is true that 
magnetically polarizable atoms can assemble to form a 
medium with nontrivial permeability. However, natural 
magnetism is a quantum effect which does not disappear 
at zero frequency. The so-called artificial magnetism is a 
conceptually different and a purely classical phenomenon. 
In particular, it should be noted that while the magnetic 
polarizability of a particle is always defined, and can, in- 
deed, be large if Ohmic losses are sufficiently small, the 
effective permeability of a medium is not always defined. 
The results of Ref. [lfl indicate that the effective parame- 
ters can only be introduced far away from magnetic res- 
onances. 

Silveirinha's correctly writes that the quadratic fluctu- 
ations of the field of the type (<5e x <5b) (e and b are the 
"microscopic" electric and magnetic fields) can not be ne- 
glected in the general case. However, in Ref. [j, I was not 
concerned with the general case but only with the case 
of electromagnetically homogeneous media which can be 
reasonably characterized by effective medium parame- 
ters. In such media, the quadratic fluctuations can be 
and should be, in fact, neglected. Accounting for such 
fluctuations amounts to exceeding the precision of the 
underlying approximation (that is, the approximation in 
which the medium is characterized by spatially-uniform 
effective parameters). Silveirinha did not demonstrate 
that the quadratic fluctuations can be large while the 
medium is still electromagnetically homogeneous accord- 
ing to the criteria stated above. And the quantitative 
arguments given in this Comment are quite indicative of 
the contrary. 

E) Unsubstantiated criticism. Silveirinha's criti- 
cism of Ref. |6| which appears after Eq. (S62) is not based 
on any scientific arguments. To be more specific, Sil- 
veirinha has derived (S62) and stated, essentially, that 
it follows from (S62) that Ref. is erroneous as shown 
by Efros. However, the argument of Ref. [H is based on 
the traditional expression for the heating rate (the one 
equivalent to (S62)) and the alternative expression is only 
mentioned as a possibility. Also, the arguments of the 
earlier paper by Efrosi^ are considered and discussed in 
Ref. at length. Finally, the main subject of Ref. are 
natural diamagnetics. It almost appears that Silveirinha 
did not look at the content of Ref. [g at all. 

In summary, I have shown that the main conclusions 
of Silveirinha's paper are based on an unphysical excita- 
tion model in which an external current which overlaps 
with the medium but is not subject to constitutive re- 
lations is used. However, if we use instead the conven- 
tional and experimentally-relevant model of excitation by 
external fields, we would find that Silveirinha's formulas 
contain the very results he wanted to disprove. I have 
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demonstrated this explicitly for the heating rate, q. I did 
not consider the density of electromagnetic energy as this 
quantity is not measurable experimentally and is, there- 
fore, irrelevant (only the total electromagnetic energy of 
a body is measurable, but this quantity is the same in 
the conventional theory and in the theory of Ref. 0) ■ As 
for the Poynting vector, S, in a steady state and, particu- 
larly, for monochromatic fields, it must satisfy the conti- 
nuity equation q + V • S = (time averaging is assumed). 
I have shown that the quantity q that follows from Sil- 



veirinha's Eq. (S60) is exactly the same as I have derived 
in Ref. 0. I have also shown in Ref. that the expres- 
sion for S which contains the cross product E x H is not 
consistent with this expression for q while the expression 
which contains E x B is. Therefore, if Silveirinha adopts 
the correct excitation model and uses it to re-compute 
the Poynting vector, he would surely find an expression 
which is quite in agreement with my previous results and 
with Silveirinha's own Eq. (S60). 
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